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The conditions of optimum orientations (lowest barrier 
and largest interaction radius) for deformed colliding nuclei 
are introduced in "cold" fusion of superheavy nuclei. Also, the 
role of (octupole and) hexadecupole deformations is studied. 
We have used the proximity potential and applied our method 
to Ca-induced reactions. 

PACS number(s): 25.60.Pj, 25.70.Gh, 27.90. +b 

Collisions between the deformed, oriented nuclei have 
been of much interest from time to time. In early 1980's 
it got trigged off, beginning with a suggestion by Greiner 
[1] that oriented 238y_|_238y collisions could lead to a very 
long lived (life-time ^ 10^^" sec) giant molecule. A num- 
ber of calculations were made [2-6] which all resulted in 
showing that the barrier is lowered due to deformations 
and orientations of colliding nuclei and that it is lowest 
for the 0°-180° orientations of two ^ssy nuclei, known 
as pole-to-pole (p-p) or nose-to-nose configuration. Note 
that 238u ^ prolate deformed nucleus, and hence the 
above result is only for prolate-prolate collisions, though 
quoted in the literature loosely for all oriented collisions. 
In fact, we show in the present paper that colliding nu- 
clei with different signs of their quadrupole deformations 
result in lowest barrier for different orientations (see Ta- 
ble 1). For example, for prolate-oblate collisions, the 
barrier is lowest for 0°-90° (equator-equator crossed, in 
short, e-c) configuration, as is envisaged very recently by 
Norenberg [7]. Also, till recently [8], the role of multi- 
poles higher than the quadrupole had not been investi- 
gated. We find that the inclusion of higher multipole 
deformations is favorable for fusion in some cases only 
i.e. the barriers are lowered only for some orientations. 
Some of these results require immediate attention and 
verification by alternative methods. 

We use here the quantum mechanical fragmentation 
theory (QMFT), extended to include the higher multi- 
pole deformations and orientations degrees of freedom. 
In QMFT, cold synthesis of new and superheavy nuclei 
was first proposed by one of us and collaborators [9-11], 
where a method was given for selecting out an optimum 
"cold" target-projectile (T-P) combination. Cold com- 
pound systems were considered to be formed for all those 
T-P combinations that lie at the bottom of the minima in 
the potential energy surface of a given compound nucleus, 
calculated for all possible T-P combinations, referred to 
as " cold reaction valleys" or reaction partners leading to 
"cold fusion" [10-14]. This information on "cold fusion 
valleys" was further optimized [11] by the requirements 
of smallest interaction barrier, largest interaction radius 
and non-necked (no saddle) nuclear shapes, identifying 



the cases of "cold", "warm/ tepid" and "hot" fusion re- 
actions. The key result behind the cold fusion reaction 
valleys is the shell closure effects of one or both the re- 
action partners. The QMFT was advanced as a unified 
approach both for heavy ion collisions, leading to fusion, 
and fission of nuclei including the cluster radioactivity 
(see e.g. the reviews in [15] and the references therein). 

We choose to apply our method to the recent ex- 
periments of highly neutron-rich '^^Ca beam bombarded 
on neutron-rich actinides ^^^Th, "^^^U, ^^^'"^^^Pu and 
^**^Cto, forming the compound systems ^^"110*, ^^^112*, 
290,292 214* aj^(i 296^15* [16]. In these reactions, for near 
the Coulomb barrier energies, the compound nucleus ex- 
citation energy E* ^30-35 MeV, in between the one for 
cold (10-20 MeV) and hot (40-50 MeV) fusion reactions. 
The use of neutron-rich (radioactive) nuclei is essential 
for overshooting the centre of island of superheavy nuclei 
(the next doubly magic nucleus) and their deformations 
and orientations could provide an added advantage since 
the fusion barrier gets lowered, or, in other words, the ex- 
citation energy of compound system gets reduced. This 
means a possibility that the "warm" and/or "hot" fu- 
sion reactions could also be reached in "cold" fusion, as 
is found to be the case here in the following calculations. 

The QMFT is worked out in terms of the mass (and 
charge) asymmetry 77=(Ai-A2)/(Ai-|-A2) (and riz={Zi- 
Z2)/(Zi+Z2)), the relative separation R, the deforma- 
tions (so far A=2 only, the quadrupole deformations) 
of two nuclei (i=l,2) or, in general, the two fragments, 
and the neck parameter e [17-20]. We introduce here 
the higher multipole deformations A=3 and 4, i.e. the 
octupole and hexadecupole deformations, as additional 
new parameters. Also, two orientation angles 9i are in- 
cluded, as in [6] (see Fig. 1, illustrated for quadrupole 
deformations). So far, the time-dependent Schrodinger 
equation in rj is solved for non-oriented collisions and for 
weakly coupled rj and r]z motions: 



i7*(77,i) =i?i^*(77,i), 



(1) 



with R(t) treated classically, and P2i and e fixed by min- 
imizing the collective potential V(R,77,ryz,/32i,e)- Eq. (1), 
solved for a number of heavy systems [19,20], shows that 
a few nucleon to a large mass transfer occurs for T-P 
combinations coming from outside the potential energy 
minima, whereas the same is zero for a T-P referring to 
potential energy minima. This means that for cold re- 
action partners, the two nuclei stick together and form a 
deformed compound system. A few nucleon transfer may, 
however, occur if a "conditional" saddle exists [21]. The 
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solution of Eq. (1) is very much computer-time consum- 
ing, and hence the following (next paragraph) simplifica- 
tions are exercised on the basis of calculated quantities. 

The potentials V(R,r7) and V(R,77x) for non-oriented 
nuclei, calculated within the Strutinsky method by using 
asymmetric two-center shell model (ATCSM) , show that 
the motions in both rj and r]z are much faster than the 
R-motion. This means that these potentials are nearly 
independent of the R-coordinate and hence R could be 
taken as a time-independent parameter. This reduces 
Eq. (1) to the stationary Schrodinger equation in 77, 



where the specific nuclear surface tension coefficient 
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+ V{r,)}^''{r,) = E'^^'^i'n). (2) 



Here R is fixed at the post saddle point, a choice justified 
by many calculations [15], and by an explicit, analytical 
solution of time-dependent Schrodinger equation in riz 
coordinate [22]. An interesting result of these calcula- 
tions is that the yields (cx| or [ ^'(7/^)1^, respec- 
tively, for mass or charge distributions) are nearly insen- 
sitive to the detailed structure of the kinetic energy term 
in the Hamiltonian which consisted of Cranking masses 
B^^ consistently calculated by using ATCSM. In other 
words, the static potential V(r/) or y{r]z) contains all the 
important information of a colliding or fissioning system. 

Since the potential V(?7,R) is nearly independent of the 
choice of R- value, for oriented nuclei, we define it as sum 
of two binding energies, and the deformation and orien- 
tation dependent Coulomb and proximity potentials: 



i=l 



(3) 



Here, the binding energies are taken from Moller et 
al. [23] for Z>8, and from experiments [24] for Z<7. The 
Coulomb and proximity potentials, with higher multi- 
pole deformations included, are obtained by following the 
works of [25] and [6] , respectively. The Coulonb potential 
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2X + 1R^+^^ ^ 



Riiai) = Roi [1 + ^ Ail^f VO] , 



(4) 



(5) 



where i?oj = 1.28^ 



1/3 



0.76 + 0.8A 
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A similar ex- 



pression is obtained by Rumin et al. [26] which differs in 
the quadrupole interaction term proportional to 
The nuclear proximity potential 



7 = 0.9517 



1 - 1.7826 ( 



MeV fm ^; the surface 



thickness &=0.99 fm; and the universal function, inde- 
pendent of the geometery of nuclear system, is 



$(so) 



i(.so - 2.54)2 - 0.0852(so - 2.54)=^ 



-3.437exp(-5%) 
respectively, for sq < 1.2511 and > 1.2511, with 
So = [R- Ri{ai)cosipi - R2{a2)cosip2]/b, 



(7) 



(8) 



Vp = 47ri?76$(so), 



(6) 



the separation distance between the colliding surfaces, 
parallel to R taken along the collision axis, in units of b. 
For So to be minimum , i.e. dso/dai = 0, it follows from 
Fig. 1 that tpi = 9i — ai, tp2 = 180 — ^2 — 0^2 andtan^j = 
—R^{ai)/Ri{ai). Finally, R, the mean curvature radius 
characterizing the gap, for nuclei lying in the same plane, 
is l/R"^ = 1/ RiiRx2 + l/ R21R22 + II R11R22 + I/ R21R12, 
where the four principal radii of curvature Rn and i?i2, at 
the points D (denoted 1) and E (denoted 2) of minimum 
So, are given by Eq. (15) in Ref. [6]. For further details, 
see [6] . Recently, Misicu and Greiner [8] have also derived 
the heavy ion interaction potential by using a multipole 
expansion of the densities in a double folding proceedure. 
Such a proceedure is shown [27] to depend strongly on 
the number of terms included in the expansion. Three 
terms are found to be sufficient for the internal region of 
the nuclear potential, whereas up to five terms are shown 
necessary for the physically more relevant surface and tail 
region for heavy ion collisions. 

For the fixed orientations, the charges in (3) are 
fixed by minimizing the potential V{R,r],r]z, /3\i,di) in 
r]z coordinate (which fixes the deformation coordinates 
P\i also). In Eq. (8), for fixed R, Sq is different for 
different orientations, and for fixed Sq, R is different for 
different orientations which is used here in the following. 

Table 1 gives the orientations of nuclei for the lowest 
barrier, for all possible combinations of different signs of 
their quadrupole deformations (prolate, oblate or spher- 
ical). These barriers also lie at the largest interaction 
radii. Fig. 2 illustrates our result for prolate-oblate 
238p-j^_l_48^j.^286-|^-|^2 reaction (see solid lines, where de- 
formations are included to all orders, A=2,3,4). Accord- 
ing to the QMFT [11], as already stated above, the above 
conditions are for an optimum cold fusion reaction. In 
other words, the orientaions in Table 1 are the optimum 
orientations for cold fusion reactions using deformed nu- 
clei. We further notice from Fig. 2, that the inclusion 
of higher multipole deformations is not always favorable 
for fusion (compare solid lines, with dotted ones for A=2 
alone): the addition of (34i term (/33j=0) lowers the bar- 
riers for some sets of oientations whereas it raises them 
for the other sets of orientations (illustrated in Fig. 2 for 
two cases each). Thus, the choice of nuclei having oc- 
tupole and hexadecupole deformations for (cold) fusion 
reactions must be made judiciously, depending on not 
only the signs of their quadrupole deformations but also 
the orientation angles. 
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Fig. 3 shows the fragmentation potentials for opti- 
mum orientations of the different T-P combinations, at 
a fixed separation so=l-5 fm, forming the same com- 
pound nucleus ^^^112*. The case of spherical nuclei 
[29] is also plotted for comparisons. Apparently, due 
to deformation and orientation degrees of freedom, all 
the potential energy minima are lowered, some new min- 
ima have appeared and some old ones have disappeared. 
Specifically, new deep minima occur at "''^Cr+^'^^Ra and 
^°^Mo-|-^®°Yb, which are in addition to the ones referring 
to the region of cluster radioacticity and/ or "hot" fusion 
(involving light nuclei of masses <30). The minima that 
have disappeared refer to well known cases of Ca and 
Pb (or neighbouring) nuclei; here 48,50Qg^^238,236u g^j^j 
^°Ge-f^°^Hg combinations. Similar results are obtained 
for the compound systems 290,292^^4* g^^^ 296 ^^g* ^or 

^*°110*, however, Ca and Hg minima are still deep and 
could be used as cold fusion reactions (details to be pub- 
lished elesewhere). We further notice from Fig. 2 that, 
w.r.t. the g.s. energy, all minima now refer to much 
smaller excitation energies, some of them lying even be- 
low it. For example, for Ca minima it is reduced from 
^35 MeV for spherical nuclei to <20 MeV for deformed 
and oriented collisions. This means that, as compared to 
spherical nuclei, oriented collisions result in cooler fusion 
reactions. Furthermore, all the optimum cold oriented 
collisions involve radioactive nuclei. 

The above results are seen better in the calculated 
mass distribution yields Y{Ai) —\ '^{ri{Ai)) p 
for ly = 0. Here, are the solutions of Eq. (3) and 
BjjTj are the classical hydrodynamical masses [30]. We 
take the view that, since fragments related to the min- 
imum in y(r?) are more probable, the yields must give 
the intermediate (two) fragment formation yields or, in 
short the formation yields for a cool compound nucleus 
[29], where the contribution of barrier penetration is not 
included. Evidently, for oriented collisions, the yields for 
new T-P pairs ^^Cr-F^soRa and lo^Mo-hisoYb are larger 
than for their neighbouring Ca and Hg induced reactions. 
The (near) symmetric combination has the largest yield, 
but they are known to form necked-in shapes, signifying 
preformation of fission fragments [11,15,29]. 

Summarizing, we have extended the QMFT for use of 
oriented collisions and inclusion of higher multipole de- 
formations, which result in the reduction of excitation 
energies of the compound system formed due to different 
T-P combinations. This means that both the "warm" 
and " hot" fusion reactions could now be reached in " cold 
fusion" also. The idea of optimum orientations for cold 
fusion reactions is introduced for the first time, which 
leads to new cold fusion reaction partners. The choice 
of nuclei with hexadecupole deformations is shown to 
depend strongly on both the signs of their quadrupole 
deformations and orientation angles. 

Table 1: The optimum orientations for "cold" fusion 
of nuclei with all possible combinations of deformations. 



Here the spherical nuclei (denoted by f) are considered 
to have zero octupole and hexadecupole deformations. 
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Figure Captions 
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Fig. 1 Schematic configuration of two axially symmetric 
deformed, oriented nuclei, lying in the same plane 
(<^ = 00). 

Fig. 2 Scattering potentials for the prolatc-oblatc 
238pu+48Ar^286]^]^2*, at different orientations. 
The R-values at the top of the three lowest lying 
barriers are also shown. 

Fig. 3Fragmcntation potentials of ^*^^112* for the op- 
timum orientations of different T-P combinations 
with A=2,3,4 (solid line with symbols) and for 
spherical nuclei (solid line). For Z<7, the /32i are 
from relativistic mean field calculations using TM2 
force [28], and for Z>7 from [23]. The /Jg^ = /34» = 
for Z<7. For the spherical case, P2i = /^Si = p4i = 
0. The g.s. is the ground state energy. 

Fig. 4Calculated yields of ^^^112 for optimum orienta- 
tions of different T-P combinations, with A=2,3,4 
(solid line with symbols) and spherical nuclei (solid 
line). 
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